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ABSTRACT 



Strong gravitational lensing by groups or clusters of galaxies provides a powerful technique to measure the dark matter properties of 
individual lens galaxies. We study in detail the mass distribution of the satellite lens galaxy in the group-scale lens SL2S J08544-0121 
by modelling simultaneously the spatially extended surface brightness distribution of the source galaxy and the lens mass distribution 
using Markov chain Monte Carlo methods. In particular, we measure the dark matter halo size of the satellite lens galaxy to be 
6.0^ 2 q kpc with a fiducial velocity dispersion of 127+ 2 i, km s _1 . This is the first time the size of an individual galaxy halo in a galaxy 
group has been measured using strong gravitational lensing without assumptions of mass following light. We verify the robustness 
of our halo size measurement using mock data resembling our lens system. Our measurement of the halo size is compatible with 
the estimated tidal radius of the satellite galaxy, suggesting that halos of galaxies in groups experience significant tidal stripping, a 
process that has been previously observed on galaxies in clusters. Our mass model of the satellite galaxy is elliptical with its major 
axis misaligned with that of the light by ~50deg. The major axis of the total matter distribution is oriented more towards the centre 
of the host halo, exhibiting the radial alignment found in N-body simulations and observational studies of satellite galaxies. This 
misalignment between mass and light poses a significant challenge to modified Newtonian dynamics. 
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1. Introduction 

Observations to date indicate that the luminous parts of galax- 
ies are embedded in dark matter halos, but properties of the 
dark matter halos are difficult to probe due to the scarcity of 
lumi nous dynamical tracers beyond the vi sible parts of galaxies 
(e.g., ISofue&Rubinll200U iBenderet al.lll994[ and references 
therein). Gravitational lensing, a phenomenon where foreground 
mass distributions deflect the light rays from background source 
galaxies, provides a powerful tool to study matter distributions 
in the Universe because lensing is sensitive to the total matter 
distribution independent of the light profiles of the foreground 
lenses. Strong gravitational lensing occurs when a single source 
galaxy is lensed into multiple images; in contrast, weak gravi- 
tational lensing corresponds to weak distortions in the shape of 
the single lensed image of the background source. 

Weak galaxy-galaxy lensing has been successfully em- 
ployed to study the dark matter halo properties of galaxies both 
in the field (e.g., [B rainerd et al. 1996; dell' Antonio & Tyson 
1996 1 ; iHudson et al.1 119981: iFischer et al.l 12000 : Smith etal. 



2001 



Hoekstra et all I2003L [2 004: Mandelbau m et al.ll2006allbF 
Parker etal. 2007; Tian et al. 2009) and in clusters (e.g., 
Nataraian & Kneib 1997; Nataraian et al. 2002; Limousin et al. 



* Based in part on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. These observa- 
tions are associated with program 10876. 



120071: iNatarai an et al.ll2009l) . In this method, the tangential shear 
signals from the distorted shapes of background galaxies around 
foreground galaxies are stacked together to infer the mass dis- 
tributions of the foreground lens galaxies. Stacking is needed 
because the shear signals around individual foreground galaxies 
are typically too weak to constrain the mass distributions of indi- 
vidual lens galaxies. Consequently, scaling relations are used to 
characterise the average properties of the foreground lens galax- 
ies in order to model the stacked shear signal. In the case of 
galaxy clusters, the cluster profile also needs to be modelled in 
detail to extract the residual shear signals around cluster galaxies 
(e.g.. iLimousin et al]|2005l) . From these studies, there is a gen- 
eral trend that galaxies in denser environments have smaller dark 
matter halo sizes. 

Strong lensing has been used to constrain th e size of dark 
matter halos in the galaxy cluster Abell 1689 (Halkol a et al.l 
2007). Scaling relations are used to relate the halo size to the 
luminosity of the galaxy members. By modelling the positions 
of the multipl e images of strongly lensed background sources, 
lHalkola et all (120071) found the halos of galaxies in Abell 1689 
to be truncated compared to galaxies of equal luminosity in the 
field. This is likely due to the tidal stripping of galaxy halos 
by the strong gravitational potential of the g alaxy cluster (e.g., 
ILimousin et alJ2009cHNataraian et alJl2009f) . 

Both the galaxy-galaxy lensing and strong lensing meth- 
ods mentioned above are statistical in that they average signals 
over many lens galaxies and study the ensemble properties of 
these lens galaxies. Recently, two studies of strong lensing clus- 
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ters have yielded halo s i ze me asurements for individual cluster 
galaxies. iRichard et alJ (2010) modelled the cluster Abell 370 
and used the multiple image positions to constrain the halo size 
of a cluster gala xy that lies in the vicinity o f a giant gravitation- 
ally lensed arc. iDonnarumma et all d2010t) analysed the cluster 
Abell 611 using the multiple image positions of lensed arcs, and 
measured the halo sizes of the brightest cluster galaxy (BCG) 
and several galaxies in the vicinity of the BCG. In both studies, 
the centroids, axis ratios and position angles of the cluster galax- 
ies have been fixed to the observed light distributions. While 
studies such as the Sloan Lens ACS Survey indicate that light 
is generally a good tracer of mass for isolated galaxy lenses in 
the field ( Koop mans et al. 2006: B olton et al. 2008a,b), there can 
be misalignments especially for galaxies in dynamical environ - 
ments such as in groups and clusters (e.g., iPereira et al.l l2008). 
If there are misalignments in the cluster galaxies in Abell 370 
and Abell 611, then fixing these lens parameter values to the ob- 
served light distribution could lead to a bias in the determination 
of the halo size. Allowing all parameters to vary, while avoiding 
the bias, would in general induce parameter degeneracies. 

In this paper, we aim to constrain the halo size of the satel- 
lite galaxy of the massive elliptical lens in SL2S J08544-0121. 
We show how we can use strong lensing to constrain individ- 
ual galaxy halo size without assumptions of scaling relations 
and without enforcement of mass following light. This tech- 
nique works on group- and cluster-scale lens systems where 
lensed arcs are comparable to, or larger than, the halo sizes of 
group/cluster galaxies. In lens systems where a group/cluster 
galaxy lies close to an arc, the morphological information of the 
arc provides powerful constraints on the mass distribution of the 
galaxy. Recently, Vegetti et al. (2010) studied the system SDSS 
J120602. 09+514229. 5, which exhibits similar lensing features 
as in SL2S J08544-0 121. They used the lensing arc to detect the 
presence of mass substructure (a visible dwarf satellite galaxy 
in this case), which demonstrates their dark matter substructure 
detection method ( Veget ti & Koopman s 2009). Properties of the 
dwarf satellite such as its mass were measured, nonetheless, the 
halo size of the dwarf galaxy was assumed to be its theoretical 
tidal radius. Our aim in studying SL2S J08544-0121 is to mea- 
sure the halo size of the satellite galaxy. 

The paper is organised as follows. We summarise the ob- 
servations of the lens system SL2S J08544-0121 in Section [2] 
and describe the reduction of the observations in Section [3] We 
present the lens modelling method and the measurement of the 
satellite halo size in Section|4] Simulations that test this method 
for the case of SL2S J08544-0121 are in Sectional Finally, we 
discuss our results in Section |6]before concluding in Section|7] 

Throughout the paper, parameter constraints are given by the 
median values with the uncertainties given by the 2.3 and the 
97.7 percentiles (corresponding to 95.4% credible intervals (CI)) 
of the marginalised probability density distributions. We assume 
a flat A-CDM cosmology with H = 70 km s _1 Mpc~' and Q A = 
1 - Qm = 0.72. From the redshifts of the lens and the source 
galaxies in Section|2l one arcsecond at the lens (source) plane in 
SL2SJ08544-0121 corresponds to 5.0 (8.5) kpc. The position 
angles of the galaxies are measured eastward from north. 

2. Observations of SL2S J08544-0121 

The lens system SL2S J08544-0121, discovered in the Canada 
France Hawaii Telescope Legacy Survey (CF HTLS), is part 
of the Strong Lensing Legac y Survey (SL2S) (Cabanac et alJ 




Fig. 1. A colour image of SL2S J08544-0121 constructed from 
the F475W, F606W and the F814W images taken with the 
HST/ ACS. The primary and satellite lenses are denoted by P and 
S, respectively. The boxes mark the multiple image systems that 
are shown in detail in Fig. [2] 



of the strong Lensing Legacy survey OLzs) (Cabanac et al. 
2007; iLimousin et alJ 12009a). De tails of the ob s ervation s of 
the lens system are described in ILimousin et a l. (2009a) and 



ILimousin et alJ (l2009bh : we summarise below the relevant ones 
for this paper. 

Hubble Space Telescope (HST) Advanced Camera for 
Surveys (ACS) images were obtained for this lens system in 
snapshot mode in three filters: F475W, F606W and F814W (ID 
10876; PI Kneib). Figure Q] shows a massive galaxy, which we 
refer to as the primary lens galaxy (labelled P), that strongly 
lenses two background source galaxies: one into a bright inner 
arc (with its corresponding images, shown in the marked boxes) 
and the other into a faint outer arc (as marked on the figure). In 
this work, we focus on the bright inner arc because the outer arc 
is too faint for (i) redshift measurements and (ii) providing use- 
ful information on the system in addition to the inner arc without 
redshift information. Discussions of the arc hereafter refer to the 
inner arc unless otherwise stated. The arc is at a distance of ~5" 
from the primary lens galaxy, making this a group-scale lens sys- 
tem. This system is special in that a galaxy, labelled S, lies near 
the arc. 

The redshifts of the primary lens and the source are Zd,p = 
0.3530 ± 0.0005 and z s = 1.2680 + 0.0003, respectively 

i I I || I 

(Limousin et al. 2009a b). The redshift of galaxy S is Zd,s - 
0.35 14 + 0.0004 (Munozet al., in preparation). The uncertainties 
in the redshift measurements are statistical only; systematic er- 
rors due to, e.g., wavelength calibration and dependence on sec- 
tion of spectrum for redshift extraction, are ~0.0005 (R. Cabanac 
2010, priv. comm.). We note that a difference of Az = 0.001 be- 
tween galaxy P and galaxy S corresponds to a line-of-sight ve- 
locity difference of ~200km s at Zd,p, which is consistent with 
the expected radial velocities in a galaxy group. Based on the 
proximity of galaxy S to galaxy P and their measured Zd,p and 
Zd,S> we interpret galaxy S as being physically associated with 
the primary lens galaxy, and we refer to galaxy S as the satellite 
lens galaxy. For the analysis, we use za-Zd,p as the lens redshift 
for the system. 

The positions of the multiple images of the background 
source show a high degree of asymmetry: the northern image 
is ~8" away from the primary lens galaxy, significantly larger 
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Table 1. Properties of the lens galaxy light distribution 



Component 


X 

(") 


y 


(") 


n 


1 


<f> 
(°) 


Primary 1 
Primary 2 
Satellite 


-0.06 
0.09 

5.27 


0.21 
-0.21 
-0.45 


2.55 
1.53 
0.41 


5.74 
2.91 

2.33 


0.95 
0.71 
0.47 


114 
137 
115 


Primary 
Satellite 


=0.00 

5.27 


=0.00 
-0.45 


2.80 
0.41 


7.97 
2.34 


0.79 
0.47 


131 
115 



Notes - The top (bottom) half gives Sersic profile parameters from 
GALFIT for a two (single) component primary lens. Columns 1 and 
2 denote the coordinate of the centroid, column 3 is the effective radius, 
column 4 is the Sersic index, column 5 is the axis ratio, and column 6 
is the position angle. 



than the ~5" of the arc. This is caused by the group environment 
of the lens system: a concentration of galaxies is located about 
~ 1 ' to the east of the system an d has similar redshifts as the lens 
system. iLimousin et al. (2009b) showed that this asymmetry can 
be exploited to study the global properties of the group. 



3. HST Image Analysis 

We obtain the standard CALACS pipeline^ reduced images 
from the ESO/ST-ECF0 Science Archive. Since single exposures 
were taken in each wavelength band, these images still contain 
cosmic rays that ar e subsequently masked using L. A. Cosmic 
(Ivan Dokkuml 120011) . The data a re further reduced using the 
MultiDrizzle package for IRAF dFruchter et al.l l2009t) to cor- 
rect for the geometric distortion. The output image pixel size 
is 0.05". 

For the reconstruction of the lensed arcs, it is important to 
remove any light contribution to the arcs from the lens galaxies. 
This is done using the GALFIT software pac kage (Peng et al. 
2002). Following, e.g., iMarshall et all d2007l) and ISuvu et al. 
(2009), the point spread function (PSF) is estimated from a star 
in the field. The regions around the lensed arcs, as well as those 
of the remaining cosmic rays, are masked for extracting the lens 
galaxy light. The light distribution of the primary lens is bimodal 
with a concentrated circular component and a smoother, more el- 
liptical component. The satellite lens and the two components of 
the primary lens are fitted simultaneously with Sersic profiles 
(Sersic ffl968l) . To investigate the light distribution of the primary 
lens as a whole, we also fit a single Sersic profile to the primary 
lens. In doing so, we mask out the core of the primary lens that 
shows clear bimodality in the light distribution. We set the origin 
of the image coordinate system to be located at the centroid of 
the single-component primary lens galaxy. The best fit values for 
the Sersic profiles are tabulated in Table [TJ The profile parame- 
ters for the satellite lens are nearly the same in both the single- 
and two-component modelling of the primary lens light. 

Despite the dwarfish appearance of the satellite galaxy when 
juxtaposed with the primary lens, the satellite galaxy is a normal 
elliptical galaxy based on its effective radius in Table [TJ a nd the 
lensing derived velocity dispersion in Secti on l6.2l (Tollerud et all 



lenstng 

Hum 



Graves et al]l2009t iGeha et al.ll200l . 



1 A pipeline for calibrating ACS data, developed at the Space 
Telescope Science Institute 

2 The European Southern Observatory (ESO), and the Space 
Telescope European Coordinating Facility (ST-ECF) 



We have chosen to use the F606W image for our analysis 
since this has the highest signal-to-noise ratio (SNR) and is least 
affected by cosmic rays near the lensed arcs. 

4. Satellite Halo Size from Lens Modelling 

In this section, we describe our lens modelling of the HST 
data for constraining the halo size of the satellite galaxy. We 
model the lens system using simply-parametrised mass profiles 
(Section l4~Tl i. and sample the posterior probability distribution of 
the lens parameters using Markov chain Monte Carlo (MCMC) 
methods (Section I4.2l i. To constrain the lens parameters, we 
use either the image positions of the multiply imaged source 
(Section 14.3b or the extended surface brightness distribution of 
the lensed source galaxy (Section [4. 4I >. This allows us to quan- 
tify the amount of additional information the extended images 
provide on the mass distribution of the lenses. 

4.1. Simply-parametrised lens profiles 

For a review on grav itational lensing, we refer the reader to 
ISchneider et al.l d2006l) . We describe the mass distribution for the 
lens system as two lens galaxies (the primary and the satellite) 
in the presence of a constant external shear. 

The primary lens galaxy is modelled as a singu- 
lar pseudoisothermal elliptic mass distribution (PIEMD; 
iKassiola & Kovnerll993l) with dimensionless surface mass den- 
sity 



*p(0i,02) = 



2$ 



where 



(l +e) 2 (1 _ e) 2 



(1) 



(2) 



(8\ , 82) are angular coordinates on the image/lens plane, and e is 
the ellipticity defined as e = (1 -q)/(l +q) with q being the axis 
ratio. The strength and axis ratio of the primary lens are bp and 
qp, respectively. The distribution is appropriately translated by 
the centroid position (#i,p, #2,p) and rotated by the position angle 
<pp. 

The satellite lens galaxy is modelled as a dual pseudo 
isothe rmal elliptical mass distribution (dPIE; Eliasdottir et al. 
|2007|) with a vanishing core radius, 



*s(0i,fc) 



2 



1 



1 



Vfl; 



+ f z 



(3) 



where bs, is the satellite lens strength, # em is defined in equa- 
tion <J3J> with gs being the axis ratio of the satellite, and t is the 
"truncation radius" of the satellite. The dPIE distribution is suit- 
ably translated by the centroid position (6*is, #2,s) and rotated by 
the position angle (p§. The radial dependence of the three dimen- 
sional mass density distribution corresponding to equation (0 



p(r) oc 



1 



r 2 (r 2 + rf) 



(4) 



where r is the three dimensional radius, r t = D&t, and D& is 
the angular diameter distance between the observer and the lens. 
Consequently, for r <K r t the mass distribution of the satellite 
galaxy is isothermal (p oc r~ 2 ), but for r » r t the mass distri- 
bution falls off as r 4 . We refer to r t as the size of the satellite 
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dark matter halo (with the corresponding angula r size as t), but 
note t hat r t is roughly the half-mass radius (e.g., Eliasdotti r et alJ 
2007). The fiducial velocity dispersion is defined as 



C"dPIE 



^s A, 
6ttZ> 



ds 



1/2 



(5) 



where c is the speed of light, and Z) s (Ai. s ) is the angular diam- 
eter distanc e between the observer (lens) and the source. The 
appendix in lEliasdottir et al.l d2007l) shows the relation between 
CdpiE and the measured velocity dispersion rx*. 

The constant external shear is parametrised by a shear 
strength y ext and a shear angle ext . In terms of polar coordi- 
nates -& and tf such that 9\ — § cos(ip) and 02 - & sin(<^>), the lens 
potential describing the constant external shear is 



1 2 

lAext(#, <fi) = r Text'? COS(2(ip - e xt)), 



(6) 



where the shear centre is arbitrary since it corresponds to an un- 
observable constant shift in the source plane. Note that /e ext = 
jV 2 i^ext is zero. The shear position angle of ext = 0° or </> ext = 
90° corresponds to a shearing along the 9\ -direction or the 02- 
direction, respectively. 

We have 13 lens parameters in total: 5 for the PIEMD (f?ip, 
#2,p> q?i <P?7 bp), 6 for the dPIE (6>i_s, 6z,s, qs, 0s, £>s, and 2 
for the external shear (yen, (p ex t)- We impose uniform priors on 
all lens parameters. However, for t, we restrict its range to be 
between 0.4" and 10". The lower limit is set by the effective 
radius of the satellite galaxy (the radius within which the inte- 
grated flux of the galaxy is half its total flux) in TableQ] Any dark 
matter contribution to the satellite galaxy would lead to a value 
for t (approximately the half-mass radius) that is larger than the 
effective radius, assuming a roughly constant mass-to-light ratio 
for the luminous matter. The upper limit of 10" is based on the 
expectation that a satellite halo size larger than the arc (~5") will 
not be constrained. 



4.2. MCMC sampling 

We follow iDunklev et all d2005h to achieve efficient MCMC 
sampling of the posterior probability density function (PDF) of 
the lens parameters and to test for convergence. The PDFs for the 
lens parameters are described in Sections [4. 3l and l4~4l Efficiency 
of the MCMC sampling hinges upon the proposal density dis- 
tribution (which encodes the information regarding the direction 
and step size for obtaining the next random step in the MCMC 
chain). Briefly, the procedure is: (i) obtain a chain with accep- 
tance rate of ~0.25 based on an initial guess of the proposal den- 
sity distribution; (ii) update the proposal density distribution by 
approximating it as a multivariate Gaussian with its covariance 
matrix set to the covariances among the parameters of the previ- 
ous chain after the burn-in phase; (iii) adjust the Gaussian widths 
of the proposal density distribution by a constant factor to obtain 
a new chain with an acceptance rate of ~0.25: a global increase 
(decrease) in the Gaussian widths of the proposal density dis- 
tribution leads to a lower (higher) acceptance rate; (iv) repeat 
(ii) and (iii) until a convergent MCMC chain is obtained. We 
find that the proposal density distribution needs updating only 
a few times to achieve efficient sampling of the posterior PDF. 
We follow the pow er-spectrum convergence test introduced in 
IDunklev et al.l (120051) . In this approach, convergence is achieved 
when the chain is drawing samples throughout the region of high 
probability (so that the correlation between successive samples 
in the chain does not bias the posterior PDF) and when the ratio 
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Fig. 2. Identification of the multiple images. The three panels (I, 
II and III) correspond to the regions marked in Fig.Q] These im- 
age cutouts are taken from the F606W image with the light from 
the lens galaxies subtracted. The triangles, circles and squares 
show the three multiple image systems. 



of the variance of the sample mean to the variance of the sample 
distribution is < 0.01 (so that statistics can be obtain with good 
accuracy). 

4.3. Lens modelling using image positions 

We identify three peaks in the source surface brightness distri- 
bution in the HST F606W image. Figure [2] shows the identified 
image positions corresponding to each of the three peaks in the 
source. The three peaks are strongly lensed into 2, 4 and 7 im- 
ages, respectively. We include the images around the satellite 
galaxy which are important fo r our study of th e satellite. These 
images were not considered bv lLimousin et a l. (2009b) since the 
focus of their study was to probe the global properties of the 
galaxy group instead of local effects due to the satellite. We esti- 
mate the uncertainties in the image positions to be 0.05", except 
for the image positions near the satellite galaxy whose uncer- 
tainties we take to be 0.15" as they are more difficult to identify. 
The image posi tions provide a tota l of 2 6 constraints. 

We follow lHalkola et al.l d2006l l2008h for the lens modelling 
based on image positions. The posterior PDF for the lens param- 
eters, collectively denoted as 77, given the image positions d pt is 



(7) 



(8) 



where JV sys is the number of multiply imaged systems (=3 in our 
case, corresponding to the 3 peaks in the source surface bright- 
ness distribution), N. is the number of multiple images in sys- 
tem j, 6ij is the observed image position, pie (//) is the modelled 



likelihood 


prior 




p(Tj\d pt ) <x P(d pt \i!) Pin) ■ 




The likelihood is 




P(d P t\i]) = =- exp 
•^pt 


, N*p N L 10. . 


-0TJ d (rj)\ 2 
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image position, o- (J - is the uncertainty in the observed image po- 
sition, and Z pt is the normalisation given by 



z Pt =(2^ 2 nrR- 



(9) 



7=1 '=1 



with 



J=l 



(10) 



Positions of the surface brightness peaks in the source plane, jS ; 
(where j = 1 . . . N sys ), are needed to predict the image positions 

$u ('?)■ ^ or a gi ven set °f values for the 1 3 lens parameters i], we 
optimise for the source positions, denoted by flj, that maximise 
equation ®. Effectively, we have marginalised over the source 
position parameters by approximating the likelihood as having a 
delta function at p.-. As described in Section |4~T1 the prior PDF, 
P{rf), is uniform with t limited to [0.4", 10"]. 

We sample the posterior PDF P(jj\d pl ) based on the MCMC 
procedure described in Section l4~2l The critical curves and caus- 
tics of the most probable lens model and the source positions are 
shown in Fig. [3] The large-scale diamond-shaped caustic curve, 
the astroid, is associated with the primary lens. The presence of 
the satellite galaxy transforms an ordinary cusp of the astroid 
into a pentagon-shaped caustic curve. Due to the singular nature 
of the satellite mass profile, there is also a "cut" (shown as a dot- 
dashed curve i n the inset), which i s the limiting case of a caus- 
tic curve (e.g.. iKormann et aljll994l) . Each of the three source 
surface brightness peaks lies in a different caustic region, and is 
thus multiply imaged different number of times. The source peak 
marked by the solid square lies outside of the main astroid caus- 
tic of the primary lens, and hence has 2 images. The solid circle 
is inside both the astroid caustic and the cut, and therefore has a 
total of 5 predicted images (the caustic crossing creates 2 addi- 
tional images and the cut crossing yields 1 more). We have ob- 
served only 4 image positions (marked by open circles in Fig. [3) 
because the 5 th predicted image, which lies near the centre of 
the satellite and is less magnified than the other images, could 
not be detected given the noise levels and the residuals of the 
lens galaxy light near the centre (see Fig.|2|. The source marked 
by the solid triangle, being inside the pentagon caustic, has two 
more predicted images than the solid circle and is thus a 7-image 
system. The creation of such a 7-image system is not possible 
without the satellite galaxy - the primary elliptical lens galaxy 
with its simple astroid can only produce 4 images (excluding the 
central, often demagnified, image that would lie near the centre 
of the primary lens). All in all, we are able to reproduce all the 
observed image positions without predicting extra images that 
can be detected in the data. The root mean square separation be- 
tween the observed and the corresponding predicted image posi- 
tions of the most probable lens model is 0.051". 

The results of the MCMC sampling for a subset of the 13 
parameters are shown in Fig. [4] There is a degeneracy between 
bp and y ext : a higher lens strength produces more internal shear 
so that less external shear is needed to explain the configuration 
of the multiple images. Similarly, there is a correlation between 
qp (not shown) and y ext : a rounder primary lens (i.e., higher qp) 
with less internal shear requires more external shear. The figure 
shows that both b§ and t are not well constrained. In particular, t 
shows a full degeneracy in the prior range of 0.4" and 10". The 
L-shaped contours of t and Z?s indicate that the posterior cannot 
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Fig. 3. Critical curves (black) and caustics (grey) of the lens 
models. Dashed (solid) curves correspond to the most probable 
lens model with image positions (extended image surface bright- 
ness) as constraints. The dot-dashed curve in the inset marks the 
cut of the satellite. The cross (plus) symbol denotes the mod- 
elled centroid position of the primary (satellite) lens. The three 
types of open symbols mark the three sets of multiple image sys- 
tems. The solid symbols are the corresponding modelled source 
positions based on the most probable lens parameters from the 
extended image modelling. The shapes of the critical and caus- 
tic curves curves are similar for the two models, but the model 
based on image positions (dashed) has a rounder primary lens 
with higher external shear. 



be well described by a multivariate Gaussian distribution. As a 
result, the multivariate Gaussian proposal density distribution is 
suboptimal and the efficiency of the MCMC sampling is ham- 
pered by these parameter degeneracies. In fact, we could not ob- 
tain a convergent chain with > 2 x 10 6 samples even after several 
updates of the proposal density distribution due to the parameter 
degeneracies. The results presented in this section are based on 
a chain of length 10 6 after the burn-in phase. 

Table [2] lists the marginalised lens parameters in the second 
column. The value of t , which nearly spans the entire range of the 
prior, shows that we cannot constrain the halo size of the satellite 
galaxy using image positions as constraints. The remedy is to use 
more information, such as the surface brightness of each pixel in 
the extended lensed images that is described next. 

4.4. Lens modelling using extended images 

In order to reconstruct the spatially extended surface brightness 
distribution of the lensed images, we need to model both the lens 
mass distribution and the extended source surface brightness dis- 
tribution. We follow Suyu et al. (2006) to reconstruct the source 
surface brightness, denoted by s, on a grid of pixels given a lens 
potential model. The reconstruction is a Bayesian regularised 
linear inversion, where we denote A as the strength of regular- 
isation and g as the form of regularisation. This reconstruction 
returns the most probable source surface brightness distribution 
Smp, and also the Bayesian evidence of source reconstruction, £, 
which is equivalent to the likelihood of the lens parameters. The 
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Fig. 4. Marginalised posterior PDF for y ext (external shear 
strength), bp (strength of primary lens in arcseconds), b$ 
(strength of satellite lens in arcseconds) and t (truncation ra- 
dius of satellite in arcseconds) based on image position mod- 
elling. The three shaded areas show the 68.3%, 95.4% and 99.7% 
credible regions. The halo size of the satellite galaxy is not con- 
strained using image positions. 

Table 2. Modelled Lens Parameter Values 
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posterior PDF for lens parameters, ;/, given the extended image 
surface brightness, d sb , is 



likelihood prior 

P(il\d sb , g) <x Pid^M '^n), 

where 

P(d sb \r],g) = & - P(d sb \i],Xg) 

= rdsP(rf sb |s,J7)P(sM,g), 



(11) 

(12) 
(13) 




Fig. 5. Marginalised posterior PDF for y ext (external shear 
strength), bp (strength of primary lens in arcseconds), bs 
(strength of satellite lens in arcseconds) and t (truncation ra- 
dius of satellite in arcseconds) based on extended image surface 
brightness modelling, plotted on the same scales as in Fig. [4] The 
three shaded regions show the 68.3%, 95.4% and 99.7% credible 
regions. For comparison, the lines show the 68.3% credible re- 
gions from Fig.|4]with image positions as constraints. In contrast 
to the full f-degeneracy from image position modelling, extended 
image modelling allows us to constrain the truncation radius of 

n S7 ;/ 

the satellite galaxy halo to be t - 1.19" „. 



and A is the opti mal regularisation determined by maximis- 
ing P(d sb \i], A, g). ISuvu et all (120061) gave the expression for 
P(rf s b|//, A, g) in their equation (19), and showed that the approx- 
imation in equation (fT2l is a valid one. 

We sample the posterior of the lens parameters with MCMC 
methods using the Bayesian evidence & of the source surface 
brightness reconstruction as the likelihood at each step of the 
chain. The priors on the satellite lens parameters are as be- 
fore (described in Section l4"Tl i. In the sampling process, we use 
30 x 30 source pixels; the typical pixel size is ~0.035". We find 
that our results are robust under changes in the number of source 
pixels, provided the source resolution is sufficient to describe the 
data (> 20 x 20). We assume the curvature for m of regularisatio n 
for the source surface brightness distribution (ISuvu et all 2006). 
which is applicable in this case given the relatively smooth sur- 
face brightness of the lensed images. Furthermore, we fix the 
regularisation constant to be A awg , the average of the optimal A 
values from a short test chain (the optimal A is different for dif- 
ferent lens parameter values). In other words, we approximate 
equation (fT2l further and use P(*/ SD |»7, A avg , g) for &. This approx- 
imation is valid because £ changes by < 1 between the optimised 
A and the fixed /t avg . 

The resulting constraints on the lens parameters are shown 
in Fig. [5] plotted on the same scales as Fig. |4] for compari- 
son. The MCMC sampling converges in <10 5 steps. The cor- 
relations among y ext , bp and qp share similar trends as in the 
results of the image position modelling, but are smaller in mag- 
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nitude. In contrast to Fig. [4] the strong degeneracy in t is broken 
with the use of the information from the extended images. The 
significantly tighter constraints based on the extended images 
are in agreement with the constraints based on image positions 
within the 68.3% credible regions. Table[2]shows the values for 
the marginalised lens parameters. In particular, the size of the 

satellite halo is constrained to be 1.19" + ' '„, corresponding to 

6.0 + Hkpc. 



r t - ".v>_ 20 . 

The most probable image and source surface brightness re- 
constructions are shown in Fig. [6] The image pixels around the 
satellite galaxy have been excluded in the modelling due to the 
imperfect lens galaxy light subtraction near the galaxy centre. 
One of the strongly lensed images is also excluded (the right- 
most triangle in panel III of Fig. |2| since it is blended with the 
residuals of the lens galaxy light subtraction. We opt to discard 
this particular lensed image to avoid including lens galaxy light 
residuals as false lensed features of the background source that 
could potentially bias our modelling. The top and the bottom-left 
panels show that the lensed features, apart from the cores of the 
lensed images, are reproduced by our model. The cores in the 
surface brightness clumps of galaxies tend to be sharply peaked; 
this leads to the visible residuals near the cores of highly magni- 
fied images due to source and image pixelisations. Nonetheless, 
this does not limit the modelling of the extended image structure 
- the reduced^ 2 is 0.75. The reconstructed source in the bottom- 
right panel indeed exhibits three peaks in its surface brightness 
distribution, thus confirming our image identification in Section 
4.31 The critical and caustic curves of the most probable lens 
parameters are illustrated in Fig. [3] 

Comparing Tables Q] and |2 we see that the modelled (pp is 
offset by ~30° from that of the observed light d istribution, as 
was previously noted by iLimousin et alJ d2009bl) . We attribute 
this to the presence of two brightness peaks in the primary lens 
galaxy; in such dynamically unrelaxed systems, it is not sur- 
prising that the observed light distribution inside the lensed arcs 
does not trace the mass probed around the arcs by strong lensing. 
N-body simulations show that dark matter halos can have the in- 
ner parts of their gravitational potential el ongated in a differen t 
direction from that of the outer parts (e.g., Havashi et al. 2007). 
Furthermore, simulations with gas indicate that the angular mo- 
mentum of the galaxy embedded in a typical dark matter halo is 
on average more aligned with the inner parts than the outer parts 
of the total mass distribution (e.g., Bett et al. 2010). In addition, 
studies of intrinsic ellipticity correlations of luminous galaxies 
suggest that the mass and the light of galaxies could be mis- 
aligne d by ~25°-35° (e.g JOtaimuraet al. 2009: iBrainerd et alJ 
2009). The modelled position angle of the mass distribution of 
the satellite is also misaligned from that of the light by ~50°, and 
we discuss this in Section ROI 

The halo size measurement is obtained by modelling the 
F606W image, selected based on its SNR and cosmic ray sub- 
traction. We could have also considered the F475W and F814W 
images. However, the F475W image has lower SNR than the 
F606W image and has several cosmic rays near the lensed arcs, 
which would make the analysis of F475W prone to systematic 
effects due to cosmic ray residuals. Compared to the F606W im- 
age, the F814W image has similarly low levels of cosmic rays 
but lower SNR. We repeat our analysis on the F814W to check 
for consistency. We find that t is constrained by the F814W data 
and agrees with the F606W value within the uncertainties; how- 
ever, the marginalised uncertainties of t from F814W is almost 
twice as large as that of F606W, as expected due to the lower 
SNR of F814W. Thus, our measurement of t from F606W is 



consistent with F814W and would not change drastically even if 
the F814W data were included. In the next section, we test the 
robustness of this measurement using simulations. 

5. Simulation Tests: Robustness of Halo Size 
Measurement 

We create 4 simulations that mimic the image morphology of 
SL2S J08544-0121 . In each case, we have an input model which 
we use to lens a background source consisting of three approxi- 
mately Gaussian peaks in its surface brightness distribution, sim- 
ilar to the reconstructed source in Fig. [6] We then convolve the 
lensed image with a Gaussian PSF of FWHM 0.08", and add 
uniform Gaussian noise. Figure [7] shows the images of the four 
simulations labelled as #1-4. The idea is to model each of the 
simulations using the same procedure as in the case of the real 
data to test the robustness of the satellite halo size measurement. 
The details of generating each simulation are as follows. 

5.1. Inputs of the simulations 

5.1.1. Simulation #1 

The input mass distribution is a PIEMD for the primary lens, a 
dPIE for the satellite galaxy, and a constant external shear. The 
parameter t for the satellite is set to 2", whereas the other pa- 
rameters of the input model are chosen so that the lensed image 
resembles SL2S J08544-0121. This simulation corresponds to 
the ideal case where the lens profiles we use for modelling are 
identical to the profiles we use to create the simulation. 

5.1.2. Simulation #2 

This simulation is similar to Simulation #1 (PIEMD + dPIE + 
constant external shear), except the input value for the truncation 
radius of the satellite galaxy is t — 5". In this simulation, the in- 
put t is approximately the extent of the lensed arc. We expect 
the constraints on t to weaken substantially as its value becomes 
comparable to, or larger than, the arc length since the morphol- 
ogy of the arc is not so sensitive to a halo size that extends be- 
yond the arc. This simulation is designed to probe the limit in 
which we can use this method to measure the halo size of the 
satellite. 

5.1.3. Simulation #3 

In this simulation, the input primary lens mass distribution 
is shallower than isothermal, as one may expect for m assive 
systems dominated by dark matter (e.g., IVegetti et alJ [2010; 
New man et all 2009). Specifically, we desc ribe the primary lens 
as a singular power-law ellipsoid (SPLE; Barkana 1998), and 
we set the slope to be y — 1.8 (such that the density p3D ~ r~ y ). 
We continue to use dPIE with t - 2" for the satellite and input a 
constant external shear. 



5.1.4. Simulation #4 



As described in Section [2] and ILimousin et alJ (I2009b). 
SL2S J08544-0121 is in a group environment and the external 
shear is mainly generated by a concentration of galaxies located 
about ~1' to the east. In this simulation, we create the simulated 
data using a singular isothermal sphere (SIS) profile for this mass 
concentration instead of a constant external shear. The input pro- 
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Fig. 6. Extended image and source surface brightness reconstruction. From left to right: top: observed image, reconstructed image; 
bottom: normalised image residual (in units of the uncertainty estimated for each pixel), reconstructed source. The morphology of 
the lensed arc is reproduced, and the reconstructed source surface brightness distribution shows three peaks. 



files for the lens galaxies are PIEMD and dPIE with t — 2". This 
simulation tests whether a constant external shear is a good ap- 
proximation for the more realistic situation where the external 
convergence also varies slightly across the lens system. 



5.2. Modelling results of the simulations 

We model each of the simulations using a PIEMD for the pri- 
mary lens galaxy, a dPIE for the satellite galaxy and a constant 
external shear, i.e., the same profiles used to model the real data 
in Section [4] We impose similar priors as for the real data on 
the lens parameters: uniform with t restricted to be between 
[0", 10"]. In each of the four simulations, the image position 
modelling provides no constraints on t, supporting the results 
in Section 14.31 For modelling the extended lensed images, we 
employ similar image masks as that of the real data (shown in 
Fig. [§}. Fig. [8] shows the marginalised PDF of Z?s and t from 



modelling the extended images in each simulation, which we 
describe below in turn. 

In Simulation #1, the top-left panel of Fig. [8] shows that the 
modelled truncation radius of the satellite is t = 1 .9 „ , which 

-0.7 

recovers the input value of 2". All the other satellite lens param- 
eters are also recovered within the uncertainties. Some of the 
other lens parameters that show strong degeneracies (e.g., y ext , 
bp, and q?) have their modelled values lie along the degeneracies 
but are offset from the input values. Nonetheless, the measure- 
ment of t is not so sensitive to these parameters since it is mostly 
degenerate with bs which is recovered within the uncertainties. 
This simulation shows that if the profiles we use to model the 
galaxies match the true mass profiles, then we can recover the 
size of the satellite halo within the uncertainties. We have tested 
this result using different image masks and source resolution - 
in all cases, the recovered t agrees with the input value within 
the errors. 
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Fig. 7. Four simulated images designed to mimic SL2S J08544-0121 but with different input mass distributions. Simulations are 
labelled as #1^1. Details of the simulations are given in Section lBTTl 



0.2 0.3 0.4 0.5 



0.2 0.3 0.4 0.5 



O 
V 

m 
u 




0.2 0.3 0.4 0.5 0.2 0.3 0.4 0.5 

b s [arcsec] 

Fig. 8. Marginalised posterior PDF of bs (strength of satellite 
lens) and t (truncation radius of satellite) for the simulations. 
The crosses mark the input parameter values. The three shaded 
areas show the 68.3%, 95.4% and 99.7% credible regions. In the 
cases where the truncation radius is constrained (without a full 
degeneracy in the prior range), the recovered value agrees with 
the input value within the uncertainties. 



In Simulation #2, the top-right panel in Fig. [8] shows that t 
is now completely degenerate with bs, and in fact with all the 



other lens parameters. This panel is based on a MCMC chain of 
length 10 6 that has not yet converged - we have not been able 
to obtain a convergent chain after > 2 x 10 6 steps, in contrast 
to the MCMC sampling of the other simulations that converges 
in <10 6 steps. Therefore, the results of this simulation suggest 
that as soon as the halo size of the satellite becomes comparable 
to the arc length, even information from the extended surface 
brightness does not allow the halo size to be determined. 

In Simulation #3, there is a mismatch in the modelled pro- 
file (PIEMD) and the input profile (SPLE) for the primary lens 
galaxy. The bottom-left panel in Fig. [8] shows that t is uncon- 
strained within the prior range. The main reason for the full 
f-degeneracy is that the model profile cannot accurately repro- 
duce the simulated arcs because the relative t hickness in the 
extended images is sensitive to the slop e (e.g., iDve & Warrenl 
I2005L lDveetal.ll2008t ISuvuet all 12010) . The image residuals 
due to the slope mismatch diminish the amount of information 
on t since these image residuals are more significant than resid- 
uals caused by changes in t. Indeed, when we model the pri- 
mary lens using a SPLE with a slope of 1.8 for comparison to 
the PIEMD model, we find that the maximum likelihood of the 
SPLE model is a factor of ~e 13 higher than that of the PIEMD 
model. Furthermore, we recover the input t (= 2") and bs values 
for the satellite (similar to the results of Simulation #1) when 
the SPLE model is used for the primary. Therefore, we also at- 
tribute the offset in the recovered satellite strength bs from the 
input value (shown in Fig. [8} to the slope mismatch. We note 
that a formal model comparison of the SPLE and PIEMD re- 
quires obtaining the Bayesian evidence of the lens models (i.e., 
marginalising over all the lens parameters), which is beyond the 
scope of this paper. This simulation shows that t is no longer 
constrained when the modelled profile is significantly different 
from the input profile in terms of the slope. This full degener- 
acy was not obtained in the case of the real data, confirming 
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that the PIEMD+dPIE+shear model is a good description for 
the SL2S J08544-0 1 2 1 system. 

In Simulation #4, the modelled value of the truncation radius 
is unconstrained within the 95.4% CI, as shown in the bottom- 
right panel in Fig. [8] however, within the 68.3% CI, it is con- 

strained to be t - 1.8" + ,,, which recovers the input value of 

-0.4 

2". Compared to Simulation #1, the sensitivity on t is weak- 
ened due to the mismatch between the constant external shear 
model and the SIS input. Since the degeneracy is not observed 
in the case of the real data, this suggests that the external con- 
vergence due to the environment varies weaker across the strong 
lens system than that of an SIS (which has ksis °c 1 /#), making 
the constant external shear a good approximation. 

The results of the simulations indicate that if the truncation 
radius is constrained, then it agrees with the input value within 
the uncertainties. 



6. Discussions 

This is the first time the halo size of a satellite galaxy is measured 
without assuming a scaling relation and without fixing any of the 
satellite lens parameters based on the observed light distribution. 



6.1. Tidal radius of the satellite 

Is the measured halo size of -6.0 kpc for the satellite galaxy rea- 
sonable? We can answer this question by comparing our halo 
size measurement to the theoretically estimated tidal radius. 

The tidal radius of the satellite is the location between the 
primary and the satellite galaxies where the gravitational tidal 
force due the primary is equal to the gravitational force of the 
satellite. Mass particles that originally belong to the satellite and 
are further from the satellite centre than the tidal radius will 
fall into the gravitational potential of the primary lens galaxy. 
By approximating the satellite and primary galaxies as spherical 
isothermal halos, we estimate the tidal radius of the satellite halo 
to be 



'"tidal 




(14) 



where r is the three-dimensional dist ance between the centre s 
of the satellite and the primary lens (Metcalf & Ma daull200lh . 
From observations and lens modelling, we have only the pro- 
jected two-dimensional distance, R, between the satellite and 
the primary. Nonetheless, we can derive the probability distri- 
bution P{r\R) by assuming that the distribution of satellites fol- 
low the total matter distribution of the primary lens halo. This 
is motivated by N-body simulations that show the distribution 
of subhal os follow roughly the mass distribution of the host 
halo (e.g., iNagai & Kravtsovl 120051; lKlvpinetal.ll2010t) . To en- 
sure that P{r\R) is normalised, we impose a hard cutoff for the 
isothermal distribution for the primary halo at r cut , i.e., we assert 
p(r) - for r > r cut . For r cut between 500 kpc and 2000 kpc (the 
expected range for a group-scale halo), we find that our PDF for 
P(r\R) is insensitive to r cut for r <k r cut . Using the derived P(r\R), 
equation 1141 . the modelled lens strengths and the projected sep- 
aration between the lenses of R — 5.1" = 25.5 kpc in Table [2] 
we show in Fig. |9]the PDF for P(r t idai|.R)- Our measurement of 
the truncation radius of 6.0^ q kpc based on the lens modelling 
in Sectionl4~4lis thus consistent with the tidal radius PDF. 
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Fig. 9. Probability density of the tidal radius of the satellite 
given the modelled projected distance of R — 25.5 kpc and 
the modelled lens strengths. Our halo size measurement of 
6.0+2Q kpc ( sna ded) is compatible with the theoretical tidal ra- 
dius. 



6.2. Galaxy halo size comparison 

We use equation (0 to derive the fiducial velocity dispersion of 
our satellite galaxy from Table [2] (t^pie = llH^fekms . Figure 
[TUl shows our results in comparison to previous studies of halo 
sizes based on strong lensing in clusters and galaxy-galaxy lens- 
ing both in the field and in clusters. We note that different stud- 
ies often employ different galaxy luminosities, scaling relations 
and reference luminosity/velocity dispersion, so the compari- 
son illustrated is more qualitative than quantitative. Our mea- 
surement, marked by an inverted triangle, shows that the halo 
size of the satellite galaxy in SL2S J08544-0121 is significantly 
smaller than isolated field galaxies of comparable velocity dis- 
persion. In fact, the satellite galaxy halo size is similar to those 
of galaxies in clusters. One may expect that galaxy members in 
group environments would have individual halos which are on 
average larger than those in cluster environments and smaller 
than isolated field galaxies. However, the primary lens galaxy in 
SL2S J08544-0121 is the brightest group galaxy (BGG), so the 
satellite galaxy is probably near the centre of the group poten- 
tial] which could make the stripping of the satellite galaxy more 
efficient. Indeed, hydro dynamical N-body simulations done by 
iLimousin et alJ d2009ch and t he ana lysis of the galaxy cluster CI 
0024+16 by Natar aian et al.1 (2009) both indicate that in galaxy 
clusters, the halo size of a galaxy is smaller the closer the galaxy 
is to the cluster centre. Our measurement of the satellite halo size 
therefore supports the tidal stripping of galaxy halos in dense en- 
vironments found in earlier lensing studies. 



3 more precisely, the satellite galaxy is near one of the possi- 
bly two centres of the group potential since the luminosity map of 
SL2SJ08544-0121 is bimodal with one peak at the BGG and an- 
other peak at a collection of galaxies that are ~Y east of the BGG 
( ILimousin et al.ll2~0 09b) 
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Fig. 10. Comparison of the galaxy halo size measurements. For 
statistical measurements, r t is the fiducial truncation radius and 
o" is the fiducial velocity dispersion. The direct measurement of 
our satellite galaxy in SL2S J08544-0121 is marked by an in- 
verted triangle, confirming the tidal stripping of galaxy halos in 
dense environments. 



6.3. Radial alignment of the satellite 



As mentioned in Section 14.41 the total mass distribution of the 
satellite is misaligned with the observed light distribution by 
~50°. Specifically, with respect to the line connecting the mod- 
elled centres of the satellite and the primary lens galaxies, the 
major axis of the modelled satellite mass distribution is at an an- 
gle of A9° +2l whereas the major axis of the light distribution is 
at 99° + 2° (95.4% CI), where both angles are measured coun- 
terclockwise. Compared to the light distribution of the satellite, 
its mass distribution is elongated more towards the centre of the 

primary, exhibiting the so-called "radial alignme nt". 

St udies based on both observat i ons (e.g., [Pe reira & Kuhn 

120051: lAgustsson & Brainerdl 1200 6: Faltenbacher et al. 2007) 
and N -body simulations (e.g., iPereira eta l. 2008; Kneb e etaO 
2008a) show that the major axes of satellite galaxies (subhalos) 
tend to point towards the centre of mass of their host. These ra- 
dial alignments are detected by averaging over many satellite 
galaxies. In the case of observations, the position angles of the 
luminous parts of the satellites are measured and compared to 
the direction of the host galaxy. In the case of N-body simula- 
tions, mass particles belonging to subhalos are identified, and 
the shapes of the subhalos are measured. By followi ng the evo- 
lution of the alignment of the subhalos with time, IPereira et al.l 
(2008) argued for tidal torquing as the mechanism for produc- 
ing the radial alignment. One would expect that dark matter ha- 
los, being more extended and having lower spins than the galax- 
ies, would be more easily torqued. Indeed, Kne be et al.l d2008b) 
found that the central parts of the subhalos were less radially 
aligned than the subhalos as a whole. Furthermore, P ereira et al.1 
(2008) found that the dark matter alignment was stronger than 
the galaxy alignment. Our measurement of the orientations of 



both the luminous and the total matter distributions of the satel- 
lite galaxy in SL2S J08544-0121 also hints at the stronger total 
matter alignment that was found for ensembles of satellites in 
earlier studies. 



6.4. Implications for modified Newtonian dynamics 

Modi fied Newtonia n dynamics (MOND) was first suggested 
by Milgrom (1983) as an alternative to the existence of dark 
matter for explaining the flat rotation curves observed in spi- 
ral galaxies. Instead of requiring additional dark matter in the 
outer parts of a galaxy, MOND postulates that Newtonian grav- 
ity breaks down in this small-acceleration region, and the ac- 
celeration experienced by a test particle is larger than the 
Newtonian value. For a review on the subject, we refer the reader 
to ISanders & McGau gh (2002). Sellwood & Kosowskvl ([200 1) 
pointed out several tests that could rule out simple theories of 
modified gravity, one of which is the misalignment between the 
axes of dark matter halos and those of the visible matter (i.e., 
mass must follow light in MOND since there is no dark matter). 
In our analysis of SL2S J08544-0121, we find that the di- 
mensionless surface mass density of the satellite galaxy, ks, 
is misaligned with that of the light by ~50°. Under General 
Relativity (GR), ks is the surface mass density of the total mass 
distribution, up to a constant factor. However, in MOND, the de- 
rived ks does not in general sca le linearly with the surface mass 
density (e.g., IZhao et all 120061) . The key question for testing 
MOND is whether the mass distribution of the satellite galaxy 
that is modelled using MOND would also be misaligned with 
the observed light distribution. iMortlock & T urner (2001) con- 
sidered gravitational lensing in MOND under the assumption 
that the deflection of photons is twice that of massive particles 
moving at the speed of light, in accordance with GR. This as- 
sumption must be valid in the Newtonian limit of any MONDian 
lensing theory. The authors find that while gravitational lensing 
deflections depend only on k in GR (thin-lens approximation), 
this is no longer true in MOND where in general the deflections 
depend on the mass distribution of the lens along the line of 
sight. Nonetheless, Mortloc k & Turner (2001) showed that for 
a typical elliptical galaxy, the thin-lens approximation holds; the 
deflection angles are the same as the GR predictions in the inner 
parts (where the projected distances a re <2kpc) and ap proach 
an asymptotic value in the outer parts. IChiu et alj ( 120 06) found 
qualitatively similar trends for the deflection angles in a rela- 
tivistic formulation of MOND by Bekenstein (2004; the tensor- 
vector-scalar theory). Therefore, modelling the satellite galaxy 
in MOND can be seen as modelling the satellite galaxy in GR 
with a different deflection angle field, or equivalently, a differ- 
ent mass radial profile. In a strongly lensed system, the position 
angle of an elliptical mass distribution is determined based on 
the observed image configuration and is quite insensitive to the 
actual radial profile of the mass distribution. Therefore, even if 
the satellite galaxy were modelled using MOND, it would still 
need to have its mass distribution elongated in a direction that is 
misaligned with the light by ~ 50° to fit the lensing observations. 
This misalignment cannot be explained by MOND and supports 
the dark matter hypothesis. 

7. Conclusions 

We modelled the mass distributions of the group-scale strong 
lens system SL2S J08544-0121 that consists of a massive el- 
liptical lens and a satellite lens galaxy. The lens galaxies are 
modelled as parametrised profiles, with the truncation radius of 
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the satellite as one parameter to characterise its halo size. Using 
Markov chain Monte Carlo methods to sample the lens mass pa- 
rameters, we compare the results based on image position mod- 
elling and based on extended surface brightness modelling. From 
this we conclude: 

- In comparison to image position modelling, extended sur- 
face brightness modelling provides significantly tighter con- 
straints on the lens parameters, although the results of the 
two approaches agree within the uncertainties. 

- Image position modelling provides no constraints on the halo 
size of the satellite galaxy. 

- Extended surface brightness modelling constrains the halo 
size of the satellite to be 6.0+^ kpc (95.4% CI). 

- Modelling of simulations that resemble the system 
SL2S J08544-0121 shows the recovery of the input trun- 
cation radius if it is constrained. Our three-component 
lens model (consisting of PIEMD, dPIE and constant 
external shear) provides an adequate description of 
SL2SJ08544-0121, 

- This technique works on galaxy group- or cluster-scale 
lenses where the galaxy member halo size is spatially less 
extended than features of the lensed arcs. 

- The measured satellite halo size is in agreement with its es- 
timated tidal radius, indicating that galaxies in group envi- 
ronments experience tidal stripping of their halos. This con- 
firms the tidal stripping of galaxy halos in dense environ- 
ments found in numerical simulations and in previous lens- 
ing analysis of galaxy clusters. 

- The position angles of the total matter distribution and the 
light distribution of the satellite are misaligned by ~50°. The 
major axis of the total matter distribution points more to- 
wards the centre of the primary halo, consistent with the de- 
tected radial alignment of satellite galaxies in previous ob- 
servational studies and numerical simulations. 

- The misalignment between mass and light of the satellite 
galaxy is a serious challenge to MOND. 

The method developed in this paper for measuring the galaxy 
halo sizes is general and applicable to other lens systems. Indeed, 
there are already several group- and cluster-scale lens systems 
known to contain strong lensing arcs with nearby group or clus- 
ter galaxy member(s). This new galactic halo ruler, when applied 
to a sample of systems, should help advance our understanding 
of galaxy evolution and of the growth and assembly of galaxy 
groups and clusters. 
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